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ABSTRACT 

The entropy distribution of the intracluster medium and the shape of its confining potential 
well completely determine the X-ray properties of a relaxed cluster of galaxies, motivating us 
to explore the origin of intracluster entropy and to describe how it develops in terms of some 
simple models. We present an analytical model for smooth accretion, including both preheating 
and radiative cooling, that links a cluster's entropy distribution to its mass accretion history 
and shows that smooth accretion overproduces the entropy observed in massive clusters by a 
factor ~ 2-3, depending on the mass accretion rate. Any inhomogeneity in the accreting gas 
reduces entropy production at the accretion shock; thus, smoothing of the gas accreting onto a 
cluster raises its entropy level. Because smooth accretion produces more entropy than hierarchical 
accretion, we suggest that some of the observed differences between clusters and groups may arise 
because preheating smooths the smaller-scale lumps of gas accreting onto groups more effectively 
than it smooths the larger-scale lumps accreting onto clusters. This effect may explain why 
entropy levels at the outskirts of groups are ~ 2-3 times larger than expected from self-similar 
scaling arguments. The details of how the density distribution of accreting gas affects the entropy 
distribution of a cluster are complex, and we suggest how to explore the relevant physics with 
numerical simulations. 



Subject headings: cosmology: theory - 
medium — X-rays: galaxies: clusters 

Introduction 



galaxies: clusters: general — galaxies: evolution — intergalactic 



Not so very long ago, people who studied clus- 
ters of galaxies were often asked how the X-ray 
emitting gas gets so hot. The answer is simple, of 
course. If radiative cooling is negligible, then grav- 
itationally driven processes will heat diffuse gas to 
the virial temperature of the potential well that 
confines it. A tougher question would have been 
to ask why the intracluster medium has the den- 
sity that it does. In order to answer that question, 
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one needs to know what produces the entropy of 
the X-ray emitting gas. 

Entropy is of fundamental importance because 
a cluster's intergalactic gas will convect until its 
isentropic surfaces coincide with the equipoten- 
tial surfaces of the dark-matter potential. Thus, 
the entropy distribution of a cluster's gas and the 
shape of the dark-matter potential well in which 
that gas sits completely determine the large-scale 
X-ray properties of a relaxed cluster of galaxies 
(see Voit et al. 2002 and references therein). The 
gas density profile p(r) and temperature profile 
T(r) in this state of convective and hydrostatic 
equilibrium are just manifestations of the under- 
lying entropy distribution. If we wish to link these 
X-ray observables to the process of cluster for- 
mation, we therefore need to understand how the 
growth of cosmic structure generates intracluster 
entropy and how processes like radiative cooling 
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and non-gravitational heating by supernovae and 
active galactic nuclei modify that entropy distri- 
bution. 

One way to approach the problem of gravi- 
tationally driven entropy generation is through 
spherically symmetric numerical models of smooth 
accretion, in which gas passes through an accre- 
tion shock as it enters the cluster (e.g., Knight 
& Ponman 1997; Tozzi & Norman 2001). If the 
incoming gas is cold, then the accretion shock is 
the sole source of intracluster entropy. If instead 
the incoming gas has been heated before passing 
through the accretion shock, then the Mach num- 
ber of the shock is smaller and the intracluster en- 
tropy reflects both the amount of preheating and 
the production of entropy at the accretion shock. 

In reality, however, the accreting gas is lumpy, 
not smooth. Incoming gas associated with accret- 
ing subhalos enters the cluster with a wide range of 
densities. There is no well-defined accretion shock 
but rather a complex network of shocks as differ- 
ent lumps of infalling gas mix with the intraclus- 
ter medium of the main halo. Yet, despite this 
complexity, numerical simulations of hierarchical 
structure formation show that the resulting en- 
tropy profile is similar to that found in the smooth 
accretion models (Borgani et al. 2001, 2002). 

This paper outlines some simple analytical 
models designed to clarify the processes that de- 
termine the entropy of intracluster gas. As has 
become customary in this field, we will refer to 

* = ^ = j-(^y /3 — a) 

as the "entropy" of the gas, while recognizing that 
the formal thermodynamic entropy per particle for 
a gas of non-interacting monatomic particles is 

s = lnK 3 / 2 + s , (2) 

where so depends only on fundamental constants 
and the mixture of particle masses. Because we ex- 
press temperatures in energy units throughout the 
paper, Boltzmann's constant is absorbed into T 
and s becomes a dimensionless quantity. The ob- 
ject of our investigation is to understand the origin 
of the intracluster entropy distribution K(M g ), 
defined so that the inverse function M g (K) is the 
mass of gas with entropy < K. 

Section 2 computes the entropy distribution 
arising from smooth spherical accretion of cold 



gas. It presents a simple analytical formulation 
for smooth accretion relating the entropy distribu- 
tion of a cluster directly to its mass accretion rate. 
The shape of the resulting entropy distribution is 
similar to that of simulated and observed clusters 
but its normalization is too large by a factor <~2- 
3, indicating that inhomogencitics in the accret- 
ing gas must be taken into account. Because of 
this difference between smooth accretion and in- 
homogeneous accretion, the normalization of the 
intracluster entropy profile reflects the lumpiness 
of the gas that accreted onto the cluster. 

Section 3 shows how preheating and radiative 
cooling change the entropy distribution produced 
by smooth accretion. We demonstrate that mod- 
est amounts of preheating raise the entropy dis- 
tribution expected from cold accretion by an ad- 
ditive term proportional to the initial entropy of 
the incoming gas. However, large amounts of pre- 
heating, comparable to the characteristic entropy 
of the halo, suppress entropy production because 
they expand the intracluster medium and reduce 
the shock velocity at the accretion front. Includ- 
ing simple corrections for preheating and radiative 
cooling yields analytical smooth-accretion mod- 
els whose entropy distributions agree well with 
the numerical models of Tozzi & Norman (2001). 
However, the effects of preheating on hierarchical 
accretion are qualitatively different because pre- 
heating smooths the gas accreting onto a cluster, 
potentially boosting its postshock entropy much 
more than the simple additive correction applied 
to smooth accretion. 

Section 4 presents evidence suggesting that ac- 
cretion of baryons onto groups was smoother than 
accretion onto clusters. We show that groups 
must have significant entropy gradients. Other- 
wise, the observed values of core entropy cannot 
be reconciled with the observed X-ray luminosity- 
temperature (L-T) relation. Polytropic models 
consistent with both the core entropy and the L- 
T relation imply that entropy levels in the outer 
parts of groups are ~2-3 times higher than ex- 
pected from simulations without cooling or pre- 
heating and from self-similar scaling of clusters. 
These findings are consistent with existing obser- 
vations of groups and suggest a transition from 
lumpy accretion to smooth accretion below a mass 
scale - lO 14 /^ 1 M Q . 

Section 5 explores how the lumpiness of accret- 
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ing gas determines the intracluster entropy dis- 
tribution. We present a naive calculation apply- 
ing simple accretion shocks to discrete accreting 
subhalos and show that this simple picture fails 
to produce enough entropy. We then consider 
what adjustments to the preshock density and 
preshock velocity would be needed to produce the 
proper amount of entropy through simple accre- 
tion shocks. However, the situation could well be 
more complicated than this because dense accret- 
ing lumps do not necessarily thermalize all their 
incoming kinetic energy within the accreting gas. 
We therefore generalize the idea of an accretion 
shock and investigate entropy production by dis- 
sipation of turbulence and relatively weak shocks 
created as accreting subhalos circulate within the 
intracluster medium. Somehow the process of hi- 
erarchical accretion in the absence of preheating 
and radiative cooling produces self-similar entropy 
profiles in groups and clusters, and we assess the 
amount of heat input needed to preserve this self- 
similarity. If this heating mode is significant, then 
it can partially offset radiative cooling in the cores 
of clusters, and we suggest how to test for this ef- 
fect using numerical simulations. 
Section 6 summarizes our findings. 

2. Smooth Accretion of Cold Gas 

Let us first consider the case of cold accreting 
gas, in which the pressure of the incoming gas is 
negligible. We will assume a spherically symmetric 
geometry, so that mass accretes in a series of con- 
centric shells, each with baryon fraction /;,, that 
initially comove with the Hubble flow. In this sim- 
ple model, a shell that initially encloses total mass 
M reaches zero velocity at the turnaround radius 
r ta and falls back through an accretion shock at 
radius r ac . We will find that the entropy distri- 
bution in this case is determined primarily by the 
rate at which matter accretes onto the cluster and 
yields an entropy distribution between K oc M g 

and K oc Mg /3 . 

2.1. Postshock Entropy 

Because the cold accreting gas is effectively 
pressureless, the equations that determine the 
postshock entropy are 



2GM£ 



P2 



1 - 



1 



= 4pi , 



(4) 

(5) 

(6) 
(7) 



where p\ is the preshock gas density and T 2 and 
P2 are postshock quantities. Equations (6) and (7) 
are restatements of the jump conditions for strong 
shocks, assuming that the postshock velocity is 
negligible in the cluster rest frame (e.g., Landau 
& Lifshitz 1959; Cavaliere, Menci, & Tozzi 1997), 
and equation (5) is strictly true only for cosmolo- 
gies with A = 0. The postshock entropy produced 
by smooth accretion of cold gas at time t is then 



3(4 Pl ) 2 /3 
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(8) 
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where fb = 0.02h~ 2 Q^ is the universal baryon 
fraction. Thus, the entropy profile arising from 
smooth accretion of cold gas depends entirely on 
the ratio r ac /n a and the accretion history M(t). 

2.2. Shock Radius 

The ratio of the shock radius to the turnaround 
radius should remain nearly constant with time 
in the case of cold accretion. One standard ap- 
proach to estimating the virial radius of a cluster 
is to assume that it is precisely equal to half the 
turnaround radius of the shell that is currently 
accreting. Setting the shock radius equal to the 
virial radius defined in this way implies £ = 0.5 by 
definition. 

More generally, one can assume that the shock 
occurs at the radius r& within which the mean 
density is A times the critical density p CT . In 
that case, r ac = (2GM/H 2 A) 1 / 3 1 where H = 
H [Sl M (l + z) 3 + (1 - flu)} 1 / 2 in a flat cosmol- 
ogy. The corresponding turnaround radius can be 
found from the equation of motion 



GM Ar 
r= -— + ~ 



(9) 



47rr ac piw a 



(3) 



In the limit of a vanishingly small cosmologi- 
cal constant A, a bound shell obeys the familiar 
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parametric solution r = r ta [(l — cos#m)/2], t = 
tvir[(^M-sin^ M )/27r],withrta = [^GMt^/ir 2 ] 1 / 3 
for a shell that collapses to the origin at time t v - lr . 
The solution for J7a = A/3#q = 0.7 is not much 
different because the A term is always small. The 
quantity Ar/3 is never greater than O.UGM/r 2 
during the trajectory of any shell that has accreted 
by the present time, which is why we neglected 
any A-dependence in equation (5). The time to 
fall from r ac < 0.5r ta to the origin is also negli- 
gible, so we obtain r ac /r ta « (tt 2 /HH 2 A) 1 / 3 . In 
this paper, we will generally set r\ — r ta /2, unless 
stated otherwise. 

The self-similar solution of Bertschingcr (1985) 
for cold accretion with Qm — 1 and /& <C 1 pro- 
vides some support for these assumptions. In that 
model, the radius of the accretion shock remains 
fixed at 0.347 times the radius of the shell that 
is currently turning around. Because M oc t 2/>3 
in the Bertschinger solution and the shell turning 
around at time t accretes at time 2t, this model 
implies £ = 1 - 0.347 • 2 8 / 9 = 0.36. 

In a self-consistent model of a shock-bounded 
intracluster medium, the ram pressure of infalling 
gas at the accretion shock must balance the ther- 
mal pressure at the outer boundary of the hydro- 
static region (i.e., P + pv 2 must be continuous 
across the accretion shock.) Thus, the value of £ 
depends on both the accretion rate of the cluster 
and its internal structure. Appendix A develops 
self-consistent solutions for the equilibrium shock 
radius in the case of a polytropic equation of state 
and a time-varying accretion history. As long as 
the accreting gas is cold, we find that £ ~ 0.35—0.6 
for the accretion histories considered in this paper. 
However, preheating of the accreting gas can drive 
£ much lower, as we will discuss in § 3.2. 

2.3. Accretion History 

Because of the constancy of £, the entropy dis- 
tribution for cold, smooth accretion directly re- 
flects the accretion history M(t). A rough es- 
timate of M(t) can be obtained from extended 
Press-Schechter theory (Bond et al. 1991; Bower 
1991; Lacey & Cole 1993). If we restrict our atten- 
tion to virialized halos that are much more massive 
than the characteristic mass, then M oc cj~ 3 /("+ 3 ) j 
where n is the power-law slope of the pertur- 
bation spectrum and u = S c (t)D(to)/ D{t) is a 
function of the critical threshold for virialization 
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Fig. 1. — Accretion histories for halos of total mass M(to) 
at the present time to- Filled squares give the logarith- 
mic mean of the maximum progenitor mass M max in the 
merger history of a halo of mass M(to) = 1O 15 /j -1 M0, 
computed from 1000 realizations of the Lacey & Cole 
(1993) algorithm. Open squares and filled triangles show 
the corresponding quantity for M(to) = 10 14 h~ 1 Mq and 
I0 13 h~ 1 Mq. Error bars indicate one standard deviation. 
The solid line shows the best fitting parabola in log M(t)- 
logt space for M(to) = I0 15 h~ 1 Mq. The long and short- 
dashed lines show best fits for M{to) = 10 14 h~ 1 M© and 
10 13 h~ 1 Mq, respectively. The dotted line indicates the 
accretion history used by Tozzi & Norman (2001) for their 
10 15 h^ 1 Mq cluster in a ACDM cosmology. 

S c (t) and the growth function D(t) (e.g., Lacey 
& Cole 1993; Voit & Donahue 1998). In a flat 
cosmology with Qm — 0.3 at t 0} the approxima- 
tion D(t) oc t - 63 is accurate to within 6% from 
O.Olto to to- For the relevant range of power- 
spectrum indices (—2 < n < —1), we therefore 
find M cx t^, with 0.9 < ( < 1.9, yielding a power- 
law entropy distribution between K sm oc M g and 
AT sm cx Mg A . These simple relations agree well 
with those found by both spherically symmet- 
ric hydrodynamical calculations (Tozzi & Norman 
2001) and three-dimensional simulations (Borgani 
ct al. 2001, 2002). However, to do a proper com- 
parision, we need a more accurate expression for 
M(t). 

Figure 1 illustrates some useful expressions for 
M(t) derived from the merger-tree algorithm of 
Lacey & Cole (1993) for a ACDM cosmology with 
Om = 0.3, Qa — 0.7, and er 8 = 0.9. We com- 
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puted 1000 realizations of merger trees leading to 
a present day halo mass M(to), assumed that M(t) 
was equal to the maximum progenitor halo mass 
M ma i at time t, and determined the best-fitting 
parabola in logM-logt space. The coefficients of 
those best fits are given in Figure 1. As expected, 
the effective value of £ is in the range 1 < ( < 2, 
except for the lower-mass halos (10 13 /i _1 M Q ) at 
late times. Figure 2 shows the logarithmic slope 
d In K/d In M g implied by these fits. Note that this 
slope generally remains between 1.0 and 1.3 but 
rises above that in the outer parts of lower-mass 
halos because of the diminishing accretion rate as 
t approaches to- 

The near linearity of the relation between K 
and M g is what sets the effective polytropic index 
of intracluster gas outside the cluster core. Be- 
cause the underlying potential is nearly isother- 
mal, the gas density scales approximately as p oc 



-1/2 



Lg leads to P oc p lati with 7 c ff 



Thus, the relation 
1.2 and 
Observations indi- 



r , implying M g 

K OC Mg 

K cx My 1 leads to 7 e ff « 1.1. 
cate that 7 e fi ~ 1.1 — 1.2 in hot clusters (Marke- 
vitch ct al. 1998, 1999; Ettori & Fabian 1999; 
DcGrandi & Molcndi 2002), supporting the idea 
that K ~ M g when gravitationally-driven pro- 
cesses dominate the production of intracluster en- 
tropy. 

2.4. Entropy Profile 

We can now compare the entropy profile pro- 
duced by smooth accretion of cold gas to entropy 
profiles computed in other ways. To simplify those 
comparisions, we will recast the entropy profiles 
in dimcnsionless form. Because the mean mat- 
ter density within r A is Ap CI , the characteris- 
tic temperature associated with overdensity A is 
T A = GM A pm p /2r A , where M A = 4irr 3 A Ap cr /3. 
These definitions lead to a characteristic entropy 
in the baryons of K A = T A /[pm p (Afbp cr ) 2 / 3 }. 
In this section, we divide all profiles by the 
cosmology-independent entropy scale K 2 oo ob- 
tained by setting A = 200. The dimension- 
less entropy is then K = K/K200, with K200 = 
(3.3 x 10 34 ergcm 2 g- 5 /3)( M200 / 10 i5^-i MQ )2/3 

for Om = 0.3. In section 5 we will find it more 
useful to work with the characteristic entropy 
and characteristic temperature T$ obtained by 
setting r A — fta/2, but comparisons involving ob- 
servations and differing cosmologies are simpler 
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Fig. 2. — Logarithmic slope d\nK /d\nM g of the entropy 
profile K sm (M g ). The profile for a 10 15 h~ x Mq halo given 
by equation (9) is nearly linear, while the profiles for lower- 
mass halos steepen near the outskirts of the cluster. This 
steepening arises because of the diminishing accretion rate 
at late times. 



with the more definite quantities associated with 
A = 200. 

Because we are idealizing the accretion history 
as a smooth increase in M(i), there is a one-to- 
one relationship between the entropy K of a gas 
shell in the final configuration, the gas mass M g 
enclosed within that shell, and the time t that the 
shell accreted. Defining a dimcnsionless accretion 
history r)(t) = M(t)/M2oo therefore allows us to 
write the dimcnsionless entropy profile as 



K sm (r]) 



100 
~3~ 



1/3 



dlnrj 


-2/3 




2/3 


dlnt_ 




to 





(10) 



with r] = f g = M g / /feM 2 oo, so that f g is the frac- 
tion of a cluster's baryons with entropy < K. 

Figure 3 compares the dimcnsionless entropy 
profile Ksm(fg) from the smooth accretion model, 
assuming £ = 0.5 for a 10 15 h^ 1 M© cluster, with 
several other entropy profiles computed in differ- 
ent ways. The filled squares connected by a dot- 
ted line show the entropy profile from a numeri- 
cally simulated cluster, the "Santa Barbara" clus- 
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K sm (f= 0.5) 
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f g = M 4 /f b M aoo 



Fig. 3. — Dimcnsionlcss entropy profiles K(f g ) showing 
the fractional amount of gas f g = M 9 / fb M200 with entropy 
< K. Solid lines show the entropy profile derived from 
equation (10) using the accretion history for a 10 15 h~ 1 Mq 
cluster from § 2.3. Filled squares and the dotted line show 
the profile from the Bryan "Santa Barbara" cluster (SB- 
Bryan). The long-dashed line depicts the unmodified pro- 
file for an NFW halo with concentration c = 8 from Voit 
ct al. (2002). The dot-dashed line with long dashes illus- 
trates the /3-model approximation to observations described 
in the text. The short-dashed line shows the profile for a 
10 15 h- 1 Mq cluster computed by Tozzi & Norman (2001), 
and the dot-dashed line with short dashes shows the result 
of equation (10) for the TN01 accretion history in Figure 1. 



ter (Frenk et al. 1999) created by an adaptive 
mesh refinement code (Norman & Bryan 1998, 
Bryan 1999). The dark-matter density profile of 
this cluster closely follows the NFW form pdm °c 
[r(l+cr/r 20 a) 2 r 1 (Navarro, Frenk, & White 1997) 
with concentration c = 8, and the gas density also 
follows this form at r > 0.1r 2 oo- Thus, we also 
show the entropy profile of gas in hydrostatic equi- 
librium in an NFW potential with c = 8 when the 
gas density is precisely proportional to the dark- 
matter density (long-dashed line; NFW-8). This 
is the prescription sometimes used to compute the 
form of the entropy profile before it is modified 
by non-gravitational processes (e.g., Bryan 2000; 
Voit & Bryan 2001; Wu & Xue 2002a,b; Voit et 
al. 2002). It underpredicts the core entropy found 
by simulations, possibly because it does not ac- 
count for energy transfer from the dark matter to 




0.5 - 



Fig. 4. — Normalized entropy profile K/f g versus f g . 
Lines have the same meanings as in Figure 3. Note that 
the entropy profiles based on simulations and simulations 
generally follow K / f g ra 1, while those based on smooth 



accretion follow K/f g 



1.5 - 2. 



the baryons (e.g., Navarro & White 1993), and 
it overpredicts the entropy in the outer regions, 
where the cluster is not in hydrostatic equilibrium. 
In order to mimic the observed entropy profiles 
of clusters, we use a /3-model density distribution 
with (3 = 2/3 and core radius of 0.1r 2 oo (Cav- 
aliere & Fusco-Femiano 1976) with a polytropic 
relation (T oc p 7eff_1 , j e s = 1.2) relating density 
to temperature. We normalize the temperature so 
that T = T200 at the core radius, which implies 
T ~ O.5T200 at r2oo- Note that this empirically 
derived entropy profile closely corresponds to the 
numerically simulated profile without any tuning 
of f3 or 7 cff . 

The nearly linear slope of the entropy profile 
from the smooth accretion model is similar to 
that of the simulated cluster, but its normaliza- 
tion is too large. Outside the core of the simu- 
lated cluster, we find K w f g , while the smooth 
accretion model with £ = 0.5 yields K 1.6/ ff . 
Dimensionless entropy in the model of Tozzi & 
Norman (2001) has an even higher normalization, 
K w 1.8/ g for a 10 15 h~ x M Q cluster, as indicated 
by the short-dashed line (see Figure 4). Much 
of the difference in normalization between their 
numerical model and the analytical model devel- 
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Fig. 5. — Dimensionless entropy profiles K(M g ) / K20Q 
for twenty-four halos from a numerical simulation with- 
out radiative cooling or feedback. The objects were drawn 
from simulation L50+ described in Bryan & Voit (2001). 
They span a mass range from 2.5 X 10 13 /i -1 Mq to 3.5 X 
10 14 h~ 1 Mq, yet they generally have nearly linear entropy 
distributions with the same normalization found in more 
massive clusters, K(r 200) ~ K20Q 1 with no significant trend 
in mass. A bold solid line shows the median profile. 

oped in this paper stems from the different ac- 
cretion history. Their accretion law converts to 
M{t) oc ^0.931 -log t/t 0) as s hown by the dotted line 
in Figure 1. Plugging this accretion law into our 
analytical model with £ = 0.45, as suggested by 
the analysis of § 3.2 and Appendix A, produces 
much better agreement at f g > 0.5, shown by the 
dot-dashed line (i^ sm -TN) in Figures 3 and 4. The 
remaining discrepancy at f g < 0.5 comes mostly 
from the preheating assumed by Tozzi & Norman 
(2001), as we will discuss in § 3.2. 

This discrepancy between the entropy gener- 
ated by smooth accretion and the entropy pro- 
duced in simulations is even larger in groups be- 
cause of their slower accretion rates. The logarith- 
mic derivative d In r\jd In t for groups is only about 
half the large-cluster value (see Figure 1). Accord- 
ing to equation (10), the dimensionless entropy 
K sm should correspondingly be ~50% larger, ap- 
proaching ~ 3i^200 at / s ~ 1. Yet, the dimension- 
less entropy profiles of simulated groups remain 
similar to those of simulated clusters. Figure 5 
shows entropy distributions for the twenty-four 
highest-mass objects in ACDM simulation L50+ 



of Bryan & Voit (2001). Except for a few outliers, 
these halos share a dimensionless entropy distri- 
bution similar to that of the more massive Santa- 
Barbara cluster, even though their masses range 
from 2.5 x 10 13 h- 1 M Q to 3.5 x 10 14 /i" 1 M . 

2.5. Departures from Smooth Accretion 

The most obvious explanation for the differ- 
ence in normalization between the smooth accre- 
tion model and those derived from simulations and 
observations is that the accreting gas in a more 
realistic model would be lumpy, not smooth. Any 
inhomogeneity in the density of gas accreting onto 
a cluster tends to reduce the mean post-accretion 
entropy, as long as the velocity of the accreting gas 
remains unchanged. The reduction occurs because 
the postshock entropy scales as v'^ c p 1 2 ^ 3 and the 
mass-weighted mean value of p\ is larger if there 
is any inhomogeneity, anisotropy, or unsteadiness 
in the accretion flow (see Appendix B). 

This sensitivity of the postshock entropy to the 
density distribution of incoming material means 
that the entropy normalization of a cluster or 
group reflects the lumpiness of the gas that ac- 
creted onto it. In Section 4 we will show that this 
effect may actually be important — the elevated 
normalization of entropy in groups suggests that 
accretion of baryons onto groups was smoother 
than accretion of baryons onto clusters. However, 
we will first continue with our analysis of smooth 
accretion in order to include preheating and cool- 
ing in the analytical model. 

3. Smooth Accretion of Preheated Gas 

Many authors have argued that preheating of 
gas that accretes onto a cluster is needed to ex- 
plain the observed slope of the L-T relation (e.g., 
Kaiser 1991; Evrard & Henry 1991; Navarro, 
Frenk, & White 1995; Cavaliere, Menci, & Tozzi 
1997; Balogh et al. 1999; Ponman et al. 1999). 
Depending on how the analysis is done, the neces- 
sary amount of preheating ranges from Tn e 2 ^ 3 ~ 
100 keV cm 2 to over 400 keV cm 2 , corresponding to 
K ~ 10 33 erg cm 2 g -5 / 3 to 4 x 10 33 erg cm 2 g~ 5 / 3 
(but see § 4.1). If the preshock entropy level K\ is 
comparable to K sm , then we can no longer assume 
that the accreting gas is pressureless. In partic- 
ular, the Mach number of the accretion shock is 
~ (Kgm/Ki) 1 / 2 , meaning that production of post- 
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shock entropy depends on the level of preheating. 

In this section we derive an entropy jump condi- 
tion that accounts for preheating and outline how 
preheating affects the position of the shock radius, 
suppressing the postshock entropy. Because radia- 
tive cooling can offset some of the effects of pre- 
heating, we also derive a simple approximation to 
treat cooling. Then, we add both preheating and 
cooling to our simple analytical entropy profiles 
and compare them with the numerically modeled 
entropy profiles of Tozzi & Norman (2001). Be- 
cause the analytical profiles closely match the nu- 
merically computed ones, we conclude that our an- 
alytical model is a good representation of smooth 
accretion. 

3.1. Entropy Jump Condition 

A jump condition for entropy production when 
preheated gas passes through a shock can be de- 
rived from the jump conditions for other quanti- 
ties. From the density jump condition, we get 

3 A 1 



100 p 



Vin = V! - V 2 = - I 1 



M 1 



Vl 



(11) 



where V\ and v 2 are the preshock and postshock 
gas velocities, respectively, in the rest frame of the 
shock, and Vi n is the velocity of incoming gas rel- 
ative to the postshock gas. The Mach number 
M = (Spivi/bPi) 1 / 2 of the shock is then deter- 
mined by 

(M 2 - l) 2 _ l6 Pl v? n 



M 2 



15 Pi 



(12) 



Noting that Ki = P\p x 5 ^ 3 and setting v ln — w ac , 
we obtain 

(M 2 I) 2 4 8 / 3 K sm 

M 2 -~5~~KT ' ( j 

Solving this quadratic equation for the larger root 
gives 

M 2_ 48/3 K sm ( 1 + q K + VI + 2q K \ 
5 K x \ 2 J 

where q K = 10ifi/4 8 / 3 if sm . 

Now we can express the jump conditions in 
terms of K sm /K\, 

Pi 



P2 4 5/ 3 ^sm (l + y/\ + 2q K , 2q K 



2 

Zqx 



+ 



<± = ± 1 + - 
92 4 V 1 + q K + VI + 2q K 




100 



Fig. 6. — Modification of postshock entropy owing to 
preheating. Accretion of cold gas produces a postshock 
entropy K sm = t)^ c /3(4pi) 2 / 3 in the high Mach-numbcr 
limit. Preheating that raises the entropy of incoming gas 
to a preshock value K\ weakens the shock but raises the 
postshock entropy K2 above K BU1 by an amount 0.84A"i . 
The solid line in the upper panel shows K2 fK\ as given by 
equation (15). The dotted line shows the postshock entropy 
K BU1 /Ki one finds if preheating is ignored. The dashed 
line shows an approximate correction for preheating: K2 ~ 
Ksm + 0.84Ki . The lower panel more precisely shows the 
difference between K2 and K Brn . 



which lead to 



,' l + Vl + 2 gJ f 2q K 
J<2 = K sm I h — 



X 1 + 



3qx 



l + q K + ^l + 2q 



K 



5/3 



(15) 



In the limit (js; C 1, equivalent to K sia /Ki ^> 
0.25, this expression reduces to K2 ~ K sm + 
OMKi (see Figure 6; Dos Santos & Dorc (2002) 
arrive at a similar approximation following a dif- 
ferent route.) From this point of view, simply 
adding a constant value to the entropy profile pro- 
duced without preheating, as in the shifted models 
of Voit et al. (2002), seems like a good represen- 
tation of the effects of preheating. 

3.2. Preheating and the Shock Radius 

If the entropy of accreting gas substantially ex- 
ceeds that which would be created by a shock, then 
we can no longer assume £ w 0.5. Because the 
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pressure just inside the shock front must balance 
the ram pressure of accreting gas, the position of 
the shock radius depends on both the current ac- 
cretion rate and the internal structure of the clus- 
ter. When extra entropy, over and above that 
produced by cold accretion, is present within the 
cluster, the equilibrium shock radius is larger. If 
the extra entropy substantially exceeds K sm , then 
the equilibrium radius becomes so large that shock 
heating is essentially turned off. 

Appendix A presents an approximate analyt- 
ical solution for the equilibrium shock radius in 
terms of the associated shock-velocity parameter 
£ that applies when preheated gas smoothly ac- 
cretes onto a cluster. Figure 7 shows the rela- 
tion between £ and M(t) arising from the accre- 
tion histories and concentration parameters used 
by Tozzi & Norman (2001). As in the numer- 
ical models of Tozzi & Norman (2001), the ini- 
tial shock radius greatly exceeds the virial radius 
{r\ = r ta /2) because the entropy of preheating, 
K i = 3 x 10 33 erg cm 2 g -5 / 3 , is much larger than 
the characteristic entropy of the nascent dark- 
matter halo. Eventually, the halo's gravity over- 
comes the effects of preheating, pulling the shock 
radius inward to a stable value of £. The more 
massive cluster (1O 15 /i _1 M ) achieves a larger 
value of £ w 0.45 primarily because its faster ac- 
cretion rate is more effective at compressing the 
gas internal to the shock front. In the less massive 
cluster (10 14 /i _1 Mq), accretion is more gradual, 
resulting in less compression and a smaller value of 
£ w 0.35. For comparison, we also show solutions 
for £ using the accretion histories of § 2.3 but oth- 
erwise identical cluster parameters. In both cases, 
faster accretion produces a higher value of £. 

The constancy of £ in the shock-dominated 
regimes of these clusters is not trivial. On one 
hand, the rise in halo concentration with time 
tends to pull the shock radius slightly inward. On 
the other hand, the decrease in accretion rate with 
time allows the shock radius to slowly increase. In 
these particular models, the two tendencies nearly 
cancel. Note also that the model we are using is 
not a good representation for £ <C 0.5 because we 
have assumed that the time at which infalling gas 
reaches the shock radius is twice the time it took 
to reach its turnaround radius. 
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Fig. 7. — Shock-radius parameter £ as a function of the 
accreted mass r)(t) = M(t)/M'20o- The solid lines show 
solutions for £ using the model developed in Appendix A 
for the cluster parameters of Tozzi & Norman (2001). The 
dotted lines show solutions found using the accretion his- 
tories of § 2.3 but otherwise identical cluster parameters. 
Lines at the upper left show clusters of mass 10 15 Mq; 
those at the lower right show clusters of mass 10 14 Mq. 
In all cases, the initial shock radius greatly exceeds the 
virial radius = rt a /2) because the entropy of preheat- 
ing (K\ = 3 X 10 33 erg cm 2 g~ 5 ' 3 ) grossly distends the in- 
traclustcr medium. After the cluster's gravity overcomes 
the entropy of preheating, the equilbrium value of £ is de- 
termined primarily by the accretion rate onto the cluster. 



3.3. Correcting for Cooling 

Radiative cooling can significantly modify the 
entropy distribution near the center of a cluster, if 
the entropy of preheating is not excessively large. 
One illustrative way to express the effects of ra- 
diative losses is with the following equation for the 
resulting change in entropy: 



dK*' 2 
dt 



3K, 
2 



3/2 



(T) 



to 



where 
K C {T) = 



2 /n e n p \ TV 2 A(T) 



(/xm p )V2 



2/3 



t; 



2/3 



(16) 



(17) 



is the entropy level at which constant-density gas 
at temperature T radiates an energy equivalent to 
its thermal energy in the time t and A(T) is the 
standard cooling function. Figure 8 shows how 
K C (T) depends on T in the models of Sutherland 
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Fig. 8. — Threshold entropy for cooling within to = 
15 Gyr as a function of temperature and metallicity. The 
quantity K C (T) is defined in equation (17). Reading down 
from the top, the lines show K C (T) for metallicities of 
log Z/Zq = 0.5, 0.0, -0.5, -1.0, -1.5, -2.0, 3.0, and zero 
metallicity, based on the cooling functions of Sutherland & 
Dopita (1993). The single solid line highlights the value of 
K C (T) for a typical cluster metallicity of log Z/Zq = —0.5. 



& Dopita (1993) for various metallicities. The 
solid line highlights K C (T) for a typical cluster 
metallicity of log Z/Zq = —0.5. 

In a simple spherical accretion model, a gas 
shell experiences no additional heat input after it 
has passed through a shock front. Gravitational 
compression can raise the temperature of the shell, 
but its entropy does not rise. Thus, radiative cool- 
ing inexorably lowers the entropy of the shell, re- 
sulting in 



K mod ( V )= Kl' 2 { n )-K^) 



,3/2, 



where 



a r K*j\T)dt 

Zt Jt(rj) 



2/3 



2/3 



(18) 



(19) 



The term K Ta d(r]) characterizes the entropy lost 
by the gas shell that accretes at time t(rj) after it 
passes through the accretion shock. In the com- 
parisons that follow, we approximate K ra< \{rj) by 
setting T equal to = GM^m p /r ta} the char- 
acteristic temperature of the dark matter halo at 
time t, which directly ties the integrand to the ac- 



cretion history. 5 We also assume a constant metal- 
licity of log Z/Z Q = -0.5. 

3.4. Comparisons with Numerical Models 

Having developed expressions for how preheat- 
ing and radiative cooling affect the postshock en- 
tropy of smooth accreting gas, we are now in a po- 
sition to compare our analytical accretion model 
with the numerical models of Tozzi & Norman 
(2001). Their standard models assumed an initial 
entropy K\ = 3 x 10 33 erg cm 2 g -5 / 3 and a con- 
stant metallicity of 0.3Z Q . The underlying cos- 
mology was ACDM with Om = 0.3, Oa = 0.7, and 
h = 0.65, so that f b = 0.158 and t = 14.5 Gyr. 

Figure 9 shows the comparison for a 10 15 h~ x M & 
cluster. In this case, the analytical model closely 
tracks the nearly linear behavior of the numeri- 
cally computed profile at f g > 0.1. Outside the 
lowest-entropy regions, the corrections owing to 
preheating and radiative cooling are relatively 
small because K\ and JsT ra d are both less than 10% 
of the cluster's characteristic entropy. However, 
preheating is important early in the accretion his- 
tory, when the characteristic entropy of the young 
halo is much smaller, and the relationship between 
preheating and cooling has important implications 
for the fate of the central gas. 

Figure 7 shows that the first <~ 5% of the 
cluster's gas accretes nearly isentropically because 
£ <C 0.5, so its entropy remains at the initial value 
of K\ — 3 x 10 33 erg cm 2 g~ 5 / 3 . After this gas 
accretes, it is subject to radiative cooling. The 
approximate cooling model of § 3.3 gives K lad rts 
3.5 x 10 33 erg cm 2 g -5 / 3 for the earliest gas to ac- 
crete, ~ 15% higher than the entropy of preheat- 
ing. (Note that the cooling threshold if c (?20o) 01 
cluster gas at time to closely corresponds to the 
integral for K rad .) The entropy of the innermost 
gas therefore drops to zero, implying a condensed 
fraction ~ 5%, which is ten times greater than the 
<~ 0.5% condensed fraction in the Tozzi & Norman 
(2001) model. Yet, if the cooling rate in the an- 
alytical model were ~ 20% smaller, or preheating 
were ~ 20% greater, the condensed fraction would 



5 Note that this definition of characteristic temperature dif- 
fers from T200 by a small amount depending on the overden- 
sity parameter A and the halo concentration; for Um = 0.3 
and halo concentration c = 8 it equals « O.9T200 at the 
present time. 
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Fig. 9. — Entropy profiles for preheated gas that has ac- 
creted into a cluster with M(to) = 10 15 h~ 1 Mq. The solid 
line shows the numerical model of Tozzi & Norman (2001). 
The long-dashed line shows the analytical model (lf mo d) 
from this paper, computed from the accretion history as- 
sumed by Tozzi & Norman (2001), including corrections for 
preheating and radiative cooling. The short-dashed line in- 
dicates the entropy profile (K BU1 ) that results if these cor- 
rections are not applied. Dotted lines show Ki, the initial 
entropy of preheating, A" c (T2oo) the cooling threshold of 
the final cluster, and the halo entropy scale i\200- The dot- 
dashed line shows if ra d, the correction term for radiative 
cooling. Because both K\ and i^ ra d are small compared 
to i\200i these corrections have only a minor effect on the 
global entropy profile outside the very center. 

be zero. This underscores the importance of ac- 
counting for entropy produced by internal feed- 
back following cooling and condensation of intra- 
cluster gas (e.g., Voit & Bryan 2001). 

Preheating has a much more significant im- 
pact on the 10 14 /i -1 Mq cluster, yet the analyti- 
cal model remains a good match to the numerical 
model. In this case, over 30% of the gas accretes 
isentropically, while £ <§C 0.5 (see Figure 7). Even 
after the accretion shock strengthens, there is still 
a significant constant offset (rs 0. 84K\) between 
K snl and the Tozzi & Norman (2001) model. Cool- 
ing diminishes the final entropy of the gas that 
accreted isentropically, and this decrease is again 
slightly larger in the analytical model than in the 
numerical model. Note also that K c (T2oo) remains 
a good approximation to -ftT ra d for the central gas. 

The characteristic entropy of a 10 13 h" 1 Mq 





2 




I i i i i I 
Mfl ") - 1 1 4 h " 1 M 

1V1 \1-q/ — 1 VJ 11 1V1 q 


, | 

- 








T_N01 




n 

\ 






K mod 


J 

// 


m 

cm 


1 .5 




- K sm 




CM 

B 






- K rad 




o 








/ / 


erg 


1 






/ / 


■* 

n 










O 








s' K 

IY 200 




0.5 






K, 








— — 7 


K (T 200 ) 









1 






0.5 



h = M / f b M 200 

Fig. 10. — Entropy profiles for preheated gas that has 
accreted into a cluster with M(to) = 10 14 h^ 1 Mq. Lines 
have the same meanings as in Figure 9. Because the en- 
tropy of preheating is relatively significant, the accretion 
radius is initially much larger than the virial radius (see 
Figure 7), and the first third of the cluster's gas accretes 
isentropically. After a strong shock develops, the entropy 
profile exceeds K Bm by 0.84i\i. Because i\rad is not 
much less than K\, radiative cooling significantly dimin- 
ishes the central entropy. 

halo is so small that 3 x 10 33 erg cm 2 g~ 5 / 3 of pre- 
heating dominates everything. The shock radius 
never approaches the virial radius, so all the gas 
accretes isentropically. Because if ra d ~ 0.5Ki for 
the earliest gas to accrete, cooling diminishes the 
central entropy by ~ 25%. Again, the analytical 
model calls for slightly more cooling than the nu- 
merical model. 

In all three cases, our approximate analytical 
model for smooth accretion seems to represent 
with reasonable accuracy all the important pro- 
cesses operating in the numerical models. The 
analytical models show that preheating can sup- 
press radiative cooling in a given mass shell if 
K sm + 0.84^1 exceeds K TaA w K C (T 200 ). How- 
ever, the cores of these objects become isentropic 
if K\ is comparable to if 200 and larger than if ra d 
(see Figures 10 and 11). 

3.5. Inhomogeneous Preheated Gas 

What effect does inhomogeneity have on accre- 
tion of preheated gas? According to § 3.1, the 
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Fig. 11. — Entropy profiles for preheated gas that has 
accreted into a cluster with M(to) = 10 13 h~ 1 Mq. Lines 
have the same meanings as in Figure 9. In this case, the en- 
tropy of preheating exceeds the characteristic entropy of the 
final cluster, so the accretion radius never approaches the 
virial radius. Accretion is therefore virtually isentropic. As 
in the 10 13 /i -1 Mq cluster, radiative cooling significantly 
lowers the central entropy. 



postshock entropy of preheated gas is 

v? 

K 2 « ^tt^tt^ + 0.84^ (20) 



3(4 Pl ) 2 /3 

As long as the shock velocities v m are comparable 
to the accretion velocity w ac , the postshock en- 
tropy will depend primarily on the prcshock den- 
sity pi and the preshock entropy K\. When we 
consider how preheating affects the case of lumpy, 
hierarchical accretion, the second term contain- 
ing Ki is relatively straightforward to apply and 
produces a correction no greater than the entropy 
of preheating. However, if preheating can signifi- 
cantly smooth the accreting gas, substantially low- 
ering the mean mass-weighted preshock density, 
then the increase in postshock entropy owing to 
the vf n /3(Api) 2 / 3 term can be considerably larger 
than K\. We will consider some observational ev- 
idence for this effect in section 4. 

Notice also that equation (20) does not ex- 
plicitly depend on the temperature of the pre- 
heated gas. We are choosing to focus on entropy 
and density, rather than on preshock temperature, 
because this way of casting the problem implic- 



itly accounts for preshock heating owing to adia- 
batic compression. In fact, focusing on preshock 
temperature can be misleading because adiabatic 
heating actually decreases the postshock entropy 
of an accreting gas blob; compression increases 
the preshock density, lowering K sm without chang- 
ing K\. The entropy gain across the shock is 
smaller because the preshock temperature is larger 
than it would have been without the compression, 
so the heat input released by thermalization of 
the incoming kinetic energy produces less entropy. 
Thus, the postshock entropy is more closely re- 
lated to the preshock density distribution than to 
the preshock temperature distribution. 

4. Evidence for Smoothed Accretion 

The results of the previous two sections imply 
that the impact of preheating on lumpy accretion 
qualitatively differs from its impact on smooth ac- 
cretion. Entropy input preceding smooth accre- 
tion suppresses entropy production at the accre- 
tion shock because it pushes the accretion radius 
beyond the virial radius (e.g., Balogh et al. 1999; 
Tozzi & Norman 2001; § 3.2). If the entropy of pre- 
heating exceeds the characteristic entropy of the 
halo, then smooth accretion becomes virtually adi- 
abatic, leading to a nearly isentropic entropy dis- 
tribution. However, dense lumps of accreting gas 
will not necessarily be prematurely halted by the 
distended intracluster medium, so it is not clear 
that entropy production by lumpy accretion will 
be diminished in quite the same way. 

On the contrary, preheating might actually en- 
hance entropy generation over what would nor- 
mally be produced in the lumpy accretion mode. 
Because preheating tends to lower the density of 
gas within small halos, it smooths the density dis- 
tribution accreting onto larger halos, increasing 
the efficiency with which accretion shocks gener- 
ate entropy. If preheating ejects a majority of the 
gas from small halos, then this effect might be 
strong enough to cause a transition from lumpy 
accretion to smooth accretion on the mass scale 
of groups. In that case, smooth accretion onto 
groups would produce an entropy normalization 
~2-3 times higher than that expected from self- 
similar scaling of clusters. Furthermore, if the 
initial level of preheating (K\) is relatively small 
compared to the characteristic entropy of the final 
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halo (-K200), then the final halo need not have a 
substantial isentropic core. 

This section examines some intriguing pieces of 
observational evidence suggesting that preheating 
alters the entropy distributions of groups primar- 
ily by smoothing the gas that accretes onto them. 
Thus, the large amounts of preheating needed to 
explain the L-T relation in the case of isentropic 
groups might be unnecessary. First, we show that 
groups with flat entropy gradients cannot simulta- 
neously be consistent with both the observed L-T 
relation and observations of core entropy. Poly- 
tropic models of groups with an effective adia- 
batic index 7 e g w 1.2 are more consistent with 
these constraints, indicating that the entropy gra- 
dients of groups are similar to those of clusters. 6 
Then, we show that these polytropic models im- 
ply entropy levels at the virial radii of groups that 
exceed the maximum value of K (7200) ~ -^200 ex- 
pected from lumpy accretion, a conclusion sup- 
ported by recent measurements of entropy at the 
outskirts of groups (Finoguenov et al. 2002). The 
excess entropy implied by observations is simi- 
lar to the amount expected from smooth accre- 
tion onto groups, perhaps indicating a transition 
from lumpy accretion to smooth accretion at a 
mass scale ~ 10 14 M Q owing to preheating in 
smaller halos. 

4.1. Entropy Gradients of Groups 

The relationship between the X-ray luminosi- 
ties and temperatures of groups and clusters has 
long been assumed to indicate an early episode of 
preheating, but the necessary amount of preheat- 
ing is a matter of some debate. Purely gravita- 
tional structure formation calls for L oc T 2 (e.g., 
Kaiser 1986), yet observations show L oc T 2 - 5 ~ 3 
(e.g., Edge & Stewart 1991). This steepening 
of the L-T slope requires some sort of non- 
gravitational modification of the core entropy 
(Kaiser 1991; Evrard & Henry 1991; Bower 1997; 
Voit ct al. 2002). 

Some models are able to fit this relation with a 
core entropy ~ 100-150 keV cm 2 (~ 1.0 — 1.5 x 
10 erg cm 

g-s/3) (e.g., Cavaliere et al. 1997, 
1998, 1999; Voit & Bryan 2001; Bialek, Evrard, 



Direct observations of intragroup entropy that appeared 
after this paper was submitted support this conclusion 
(Mushotzky ct al. 2003; Pratt & Arnaud 2003). 



& Mohr 2001; dos Santos & Dore 2002). Other 
models require ~ 200-400 keV cm 2 (~ 2.0-4.0 x 
10 33 ergcm 3 g" 5 / 3 ) (e.g., Balogh et al. 1999; Tozzi 
& Norman 2001; Babul et al. 2002). The pri- 
mary difference between these two classes of mod- 
els is that those requiring large amounts of pre- 
heating have isentropic entropy profiles extending 
to a large fraction of the virial radius, while those 
requiring less central entropy assume a significant 
entropy gradient throughout the group. 

We can quantify this difference using the 
polytropic models described by equations (A3) 
through (A9) of Appendix A. These models re- 
quire two boundary conditions and an effective 
polytropic index 7 e ff. For this application, we set 
one boundary condition so that the pressure at 
r 2 oo is the same for all models, and we choose 
that pressure to be equal to the 7 e g = 1.2 case 
with a total gas mass /&M200 within r2oo- For the 
other boundary condition we fix the entropy at at 
0.1r2oo- Each model is therefore determined by the 
halo mass M 2 oo, the halo concentration parameter 
c 2 oo, the polytropic index j e s, and the core en- 
tropy Ko.i = if (0.1r 2 oo)- For each halo mass, we 
consider a set of polytropic models with 7 ff = 1.2 
and a set of isentropic models (7 e fj = 5/3). 

Figure 12 shows the luminosity-temperature 
relation we need to reproduce. We plot the 
quantity Lx/Tj 3 m , where 71 um is the luminosity- 
weighted temperature, because it is nearly con- 
stant at - 1.5 x 10 42 h~ 3 ergs" 1 keV~ 3 for T lum > 
1 keV and a little bit smaller below that tem- 
perature. Figure 13 shows the same quantity 
drawn from the polytropic models for two dif- 
ferent halo masses, M 2 oo = 5 x 10 13 h^ 1 M & 
and M 2 oo = 1 x 10 13 h^ 1 M Q with luminosity- 
weighted temperatures of Ti U m ~ 1.4 keV and 
Tlum ~ 0.5 keV, respectively. The halo concen- 
tration is C200 = 10 in all cases and we truncate 
the integrations for Lx and Ti um at the virial 
radius. 

Comparing Figure 12 with Figure 13 shows that 
both the isentropic models and the 7 c g- = 1.2 mod- 
els can reproduce the L-T relation, but the isen- 
tropic models call for higher levels of core entropy. 
The reason for this behavior is that the density 
profile of the isentropic models is shallower than 
p oc r~ 3 / 2 , so that Lx is dominated by emission 
from near the virial radius. The overall luminosity 
is then set by the entropy level in the outer regions 
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Fig. 12. — Luminosity-temperature relation for groups 
and clusters in terms of Lx/Tj 3 m . Solid triangles show 
measurements of clusters with insignificant cooling flows 
compiled by Arnaud & Evrard (1999). Open squares 
show cooling-flow corrected measurements by Markevitch 
(1998). Solid circles show group data from Hclsdon & Pon- 
man (2000). Note that the quantity Lx/T^ um remains 
roughly constant at ~ 1.5 X 10 42 h -3 erg s — 1 keV - 3 (up- 
per dotted line), dipping to slightly lower levels (~ 0.5 X 
10 42 /i~ 3 ergs" 1 kcV~ 3 (lower dotted line) below ~ 1 keV. 



of the group. Conversely, when 7 e g w 1.2, the den- 
sity profile is steeper than p oc r -3 / 2 . The inner 
regions then dominate the total luminosity, and 
the core entropy can be lower because a smaller 
amount of gas produces most of the emission. Be- 
cause observations of core entropy in ~1 keV halos 
indicate that ifo.i ~ 100-150 keV cm 2 (Ponman et 
al. 1999), models with significant entropy gradi- 
ents, consistent with 7 e g ~ 1.2, are probably closer 
to the actual entropy profiles of groups. 

4.2. Entropy at the Outskirts 

The preceding analysis suggests that groups are 
not isentropic but instead have entropy gradients 
more like those observed in clusters. This find- 
ing, when coupled with the elevated levels of core 
entropy observed in groups (e.g. Ponman et al. 
1999), implies that entropy at the outskirts of 
groups must also be elevated above the predic- 
tions of pure gravitational structure formation. If 
the shape of the entropy gradient is unchanged, 
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Fig. 13. — Dependence of Lx /Tj s UIn on the entropy scale 
Kq.i at 0.1r200 hi polytropic models. Tracks are shown for 
two different halo masses: 5 X 10 13 h~ 1 Mq (solid lines) 
and 1 X 10 13 h^ 1 Mq (dashed lines). Each pair of tracks 
corresponds to different effective polytropic index. The up- 
per tracks represent isentropic models (7 e ff = 5/3), and the 
lower tracks represent polytropic models with an entropy 
gradient similar to those observed in clusters (7 e g = 1.2). 
The halo concentration is C200 = 10 in all models, and 
the pressure at r2oo is the same in all models. Dotted 
lines indicate the same values of L/T 3 as in Figure 12. 
One can reproduce the observed L/T 3 ratio with cither 
7eff ~ 1-2 and A"o.i ~ 130keVcm 2 or with 7 e ff R. 5/3 and 
K .i » 350kcVcm 2 . 

then the main difference between groups and clus- 
ters must be in the normalization of the entropy 
gradient rather than its slope. 

Figure 14 shows the scaled entropy K/K 2 oo 
at r2oo implied by polytropic models having dif- 
ferent halo masses but a fixed value of core en- 
tropy _o.i = 150keVcm 2 . In polytropic mod- 
els with 7 e ff = 1.2 and this value of core entropy, 
the implied entropy at r2oo is ~ 2-4 times larger 
than if 200 over the mass range lO 13 ft- _1 M to 
10 14 h -1 Mq. For comparison, we also show the 
implied entropy for an effective polytropic index of 
7efi = 1.3, which agrees less well with the observa- 
tional constraints. In that case, the outer entropy 
levels drop to ~ 1-3 times if 200- However, the 
entropy levels at if 200 in groups simulated with- 
out preheating or cooling are rarely much greater 
than if 200, with no significant trend in mass (see 
Figure 5). Notice that these entropy enhance- 
ments are quite difficult to produce by heat input 
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Fig. 14. — Entropy at r200 in polytropic models with 
core entropy Ko.i = 150kcVcm 2 . The polytropic indices 
(7eff ~ 1-2 — 1.3) consistent with both the L-T relation and 
the observed core entropy of groups ~ 100-150 kcVcm 2 , 
imply X(r20o) ~ 2-3if200- 



alone. The value of if 200 for a 10 14 /i _1 M© halo 
is ~ 700 keV cm 2 , so doubling that entropy would 
require an enormous amount of heat. 

Direct observations of entropy at the outskirts 
of a few groups support the idea that groups have 
significant entropy gradients and higher maximum 
entropy levels than predicted by self-similar mod- 
els that do not include non-gravitational processes. 
Finoguenov et al. (2002) determined Tn e 2 ^ 3 as a 
function of radius in several groups using ASCA 
and ROSAT data and compared them with the 
entropy gradients of more massive systems. After 
dividing these gradients by M 5 g to remove the de- 
pendence on mass scale, they found that the outer 
entropy levels of groups exceeded those in clusters 
by a similar factor of 2-3, with a large scatter. 

Finoguenov ct al. (2002) interpreted this excess 
entropy as a ~ 400 keV cm 2 "entropy ceiling" pro- 
duced entirely by non-gravitational heating owing 
to galactic winds. In this interpretation, the en- 
tropy gradient results from a gradual rise in the 
entropy of the intergalactic medium external to a 
group before it accretes onto the group. For this 
level of entropy to be produced by supernovae, the 
energy injection must occur with near 100% effi- 



ciency at redshift z <~ 3 into regions of relatively 
low overdensity (~ 5p CT ). 

4.3. A Smooth- Accretion Interpretation 

We would like to offer another interpretation 
of this apparent entropy excess that is based on 
the smooth accretion models developed in this 
paper. Relatively small amounts of preheating 
(< lOOkeVcm 2 ) can eject much of the gas from 
the smaller halos that accrete onto a group or a 
large elliptical galaxy at late times, making the ac- 
creting gas more homogeneous. This more mod- 
est level of preheating can therefore break self- 
similarity by allowing the group to generate most 
of its entropy through smooth accretion. The gen- 
eral shape of the entropy profile remains similar to 
the lumpy-accretion case because it is determined 
by the accretion history (§ 2.3), but its normal- 
ization is boosted because the mean density of in- 
coming gas is reduced (§ 2.5). 

In the limit of perfectly smooth accretion, equa- 
tion (10) implies that 



K(r 2 ao) 



K- 



2.6 



200 



dlnM 
dixit 



-2/3 



(21) 



assuming that H t , 2£, and r\ are all ~ 1. The 
enhancement of entropy at r 2 oo relative to expec- 
tations from hierarchical accretion is larger for 
smaller halos because it depends on the mass ac- 
cretion rate. Inserting the d In M/d In t values at 
t = to from Figure 1 into equation (21) yields 
K(r 2 oo)/K 20 a « 3.5 for 10 13 /i^M©, w 2.4 for 
10 14 /i^M©, and w 1.7 for 10 15 h- 1 M . Note 
also that smooth accretion can increase the ef- 
fectiveness of preheating at limiting condensation 
in groups, as the initial entropy input needed to 
smooth the accreting gas is strongly amplified at 
the accretion shock. 

5. The Puzzle of Entropy Normalization 

Groups and clusters appear to have similar en- 
tropy profiles, consistent with an effective poly- 
tropic index 7 e g- ~ 1.2 at large radii. However, the 
normalizations of those profiles cannot scale with 
K 2 oo oc AfjQQ, as they would if they were deter- 
mined by purely gravitational processes. Instead, 
the relationship between entropy normalization 
and halo mass is somewhat shallower, quite pos- 
sibly because preheating has partially smoothed 
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the intergalactic medium. Thus, we would like to 
understand more quantitatively how entropy pro- 
duction depends on the homogeneity of accreting 
gas. 

The normalization of intracluster entropy in the 
smooth-accretion limit is relatively simple to cal- 
culate (see § 2), but the physics that sets the nor- 
malization during lumpy accretion is more com- 
plex. This section outlines some of the physical 
processes that might govern entropy production by 
lumpy accretion. We first present a naive model 
for merger shocks in which all the accreting gas 
belongs to virializcd subhalos and all the incom- 
ing kinetic energy is thermalized through merger 
shocks within those subhalos. This naive model 
fails to produce the required amount of entropy 
implying that either (1) the preshock density dis- 
tribution and perhaps the accretion velocity as- 
sumed in the naive model poorly represent real- 
ity or (2) some of the incoming kinetic energy is 
thermalized within the intracluster medium of the 
main halo. One way to account for the first pos- 
sibility is to adjust the preshock density distribu- 
tion, and we present an illustrative example in- 
volving power-law density distributions. The sec- 
ond possibility implies that lumpy accretion pro- 
vides a distributed source of heating throughout a 
cluster, and we conclude the section by identifying 
the hallmarks of distributed heating and suggest- 
ing how to quantify it in numerical simulations. 

5.1. Naive Merger Shocks 

The accretion process usually envisioned in 
semi-analytic models of cluster formation calls 
for all the accreting gas to lie within a subhalo 
of some kind, even if many of those subhalos sit 
below the practical resolution limit of the calcu- 
lation (e.g., Wu, Fabian, & Nulsen 2000; Bower 
et al. 2001). In a self-consistent semi-analytical 
model one would like to track how entropy devel- 
ops as these subhalos collide and merge with the 
main halo and to calculate the density configura- 
tion of the new halo from that evolving entropy 
distribution (e.g., Balogh et al. 2003, in prepara- 
tion). Instead, the density distribution is typically 
assumed to settle into a polytropic or /3-model 
configuration similar to observed density profiles, 
with a core temperature equal to the virial tem- 
perature. Here we show that accretion shocks 
within discrete, accreting subhalos do not gener- 



ate enough entropy to reproduce the characteristic 
entropy profile of a cluster. 

If the mean gas density within each of the ac- 
creting subhalos is A/f,p cr and strong shocks ther- 
malize all the kinetic energy of accretion within 
those subhalos as they cross the virial radius of 
the main halo, then for £ 0.5 we find that the 
mean postshock entropy of gas accreting at time t 



is 



K(t) 



Pi 
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K sm (t) 



(22) 
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The cosmologies and accretion histories considered 
here give K(t) ~ 0.1(A/200)- 1 / 3 K sm (t), several 
times smaller than the value found in the sim- 
ulations and observations. According to § 3.1, 
about 84% of the preshock entropy in the subha- 
los should be added to K, but this does not come 
close to fixing the problem unless the accreting ha- 
los are similar in size to the main halo, which is 
rare during the formation of a large cluster. 

This postshock entropy deficit can be recast in 
terms of the shock velocity needed to rectify it. 
Suppose that the accreting subhalos have gas den- 
sity distributions similar to that of the main halo. 
Because the gas temperature in the main halo is 
w GM [im p /2r l \, boosting the entropy distribu- 
tion of the subhalo's gas to match that of the main 
halo requires a shock velocity 



3 • 4 2 / 3 GM 
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(23) 



where w ac is the infall velocity at radius r& = 
r ta /2. A gas blob entering a cluster with an NFW 
potential would need to fall unimpeded from its 
turnaround radius to < 0.25rA before it reached 
a velocity w 2v ac . However, in the naive model 
we have outlined, infalling blobs do not fall so far 
toward the center of the main halo before being 
shocked and disrupted. If a dense lump of accret- 
ing gas does fall this far, then dissipation of its 
kinetic energy through shocks and turbulence is 
also quite likely to heat the main halo's gas. We 
consider how this mechanism affects entropy gen- 
eration in § 5.3. 
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5.2. Adjusted Merger Shocks 

If one retains the assumption that merger 
shocks thermalize all the incoming energy within 
the accreting gas, then then one must adjust the 
incoming density distribution and perhaps raise 
the shock velocity. Striking the proper blend of 
smooth accretion and lumpy accretion would then 
yield the correct entropy level. To illustrate this 
idea, we briefly outline a case in which the accret- 
ing gas fills the accreting volume instead of being 
confined to regions with mean density <~ A/;,p cr . 

Suppose that the preshock density of accreting 
gas is distributed like a power law. Define f(p)dp 
to be the fraction of accreting gas with density be- 
tween p and p + dp, so that f(p) = (p— l)/^7n/°~ p 
for a power-law index p > 1. If we then assume 
that all the gas accreting at a given time shocks 
at the same accretion radius, we can integrate 
p~ 2 / 3 f(p)dp to find the mean entropy of gas ac- 
creted at time t: 

K(t) = R ^(-^-) 2/3 KUt) ■ (24) 

V 3 \ Pmin / 

If we further assume that the accreting gas fills the 
available volume, then pi/pmi n = p/(p— 1) and 

K(t) = pV3( - P ~p 1/3 K sm (t) . (25) 

P 3 

The power-law index p is directly related to the 
density distribution within an accreting subhalo. 
In particular, for p(r) cx r~ q , we have p = 3/q. 
The case of extended singular isothermal spheres 
accreting onto the main halo therefore corresponds 
top = 3/2, and we get K w 0.89K sm {t). In the 
limit of small q and large p, the incoming density 
distribution is nearly uniform, and K w K sm (t). 
However, K grows arbitrarily small in the limit 
of q — > 3 and p — > 1, in which case the majority 
of the gas becomes concentrated at high density. 
Thus, for some value of q < 3, the mean postshock 
entropy has the desired value. 

This calculation implies that one can specify a 
preshock density distribution that yields the cor- 
rect postshock entropy distribution after passing 
through an accretion shock. However, this model 
for entropy production is not necessarily the solu- 
tion to the entropy normalization puzzle. In order 
to test it, one would need to measure the density 



distribution and infall velocity of baryonic mate- 
rial accreting onto a numerically simulated cluster 
and then compute the postshock entropy distribu- 
tion implied by this model. 

5.3. Distributed Heating 

An alternative solution to the entropy normal- 
ization problem depends on an important qual- 
itative difference between smooth accretion and 
lumpy accretion. Smooth accretion thermalizes 
all the incoming kinetic energy at the accretion 
shock, within gas that has just accreted. After 
a gas shell has become part of the intracluster 
medium, its entropy is no longer affected by the 
smooth-accretion process. However, dense lumps 
of accreting gas can penetrate the virial radius 
without substantially decelerating. They subse- 
quently move through the main cluster, stirring 
up turbulence and stimulating shocks within the 
main cluster's gas. These processes continue until 
either drag saps all of the dense lump's kinetic en- 
ergy or hydrodynamic instabilities tear the dense 
lump apart. 

The recently observed presence of cold fronts 
in clusters (Markevitch et al. 2000; Vikhlinin, 
Markevitch, & Murray 2001) suggests that dense 
subhalos can remain coherent for quite a long time 
before blending with the rest of the intracluster 
medium. Thus, much of the incoming kinetic en- 
ergy of lumpy accretion may ultimately be de- 
posited into the pre-existing intracluster medium 
as the disturbances stirred up by the accreting 
lumps dissipate and thermalize (see e.g., Gomez 
et al. 2002). If that is the case, then accretion 
provides a distributed source of heating that adds 
entropy to the entire intracluster medium, not just 
to gas that has recently accreted. 

Here we examine some crucial differences be- 
tween entropy production through smooth accre- 
tion and entropy production through distributed 
heating. First we assess the amount of entropy 
growth that lumpy accretion needs to generate in 
order to maintain the self-similar entropy profiles 
observed in simulations and in hot clusters. Then 
we examine the implications of heating the exist- 
ing intracluster medium through lumpy accretion, 
showing that the source of condensing low-entropy 
gas in the lumpy-accretion mode is different from 
its origin in the smooth-accretion mode. If dis- 
tributed heating is important, then much of the 
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core gas in present-day clusters may come from 
the cores of subhalos that accrete at ~ 0.5io- One 
can assess the importance of distributed heating 
by tracking the K(t) trajectories of Lagrangian 
gas parcels in numerical simulations after they ac- 
crete onto a cluster. 

5.3.1. Self- Similar Entropy Growth 

One notable feature of numerically modelled 
clusters in simulated CDM-like cosmologies is 
their near self-similarity. The dark-matter poten- 
tial wells approximately follow the NFW density 
profile, and the main systematic devations from 
self-similarity can be characterized with a halo 
concentration parameter that increases with de- 
creasing halo mass (e.g., Navarro et al. 1997). 
This dependence on halo mass comes about be- 
cause smaller halos tend to accumulate their mass 
earlier in time, when the mean density of the uni- 
verse is larger (Navarro et al. 1997; Bullock et al. 
2001; Eke, Navarro, & Steinmetz 2001). Numer- 
ical models including hydrodynamics also show 
that the baryonic density and temperature pro- 
files in the absence of non-gravitational heating 
and cooling processes closely track the analogous 
dark-matter profiles, except in the very center 
(e.g., Navarro, Frenk, & White 1995; Eke et al. 
1998; Frenk et al. 1999). Therefore, whatever pro- 
cesses are responsible for entropy growth during 
hierarchical structure formation act to preserve 
this near self-similarity. 

If clusters were precisely self-similar, then their 
properties would be entirely determined by their 
mass M and the time of observation, which de- 
termines the mean mass density A/? cr within the 
cluster. According to the simple spherical infall 
model of § 2.2, the turnaround radius of matter 
accreting at time t is r ta = (2GMt 2 /tt 2 ) 1 ^ 3 and it 
falls through the virial radius (r ta /2) with a kinetic 
energy per particle = GMfim p /r ta . (These 
expressions are strictly true only if A = 0, but 
they are accurate to within 10% for ACDM with 
Ha = 0.7.) Because the mean density within the 
virial radius is 3ir/Gt 2 , the entropy scale 7 associ- 



7 This entropy scale exceeds K200 by a factor (200/ A) 2 / 3 
or « 1.6 for a ACDM cosmology at the present time. We 
are using in the self-similar analysis because it depends 
explicitly on M and t alone. 



ated with a cluster's baryons is 

Kiff <Mi)2/3 



(26) 



The dimensionless entropy distribution function 
K(M g )/K c f > would then be the same for all clus- 
ters. 

One can integrate over this distribution to find 
the total classical thermodynamic entropy of a 
cluster's baryons, 



S(M, t) 



fbM 
firrip 



[ln(Mi) + const. 



(27) 



where the details of the self-similar distribution 
function are absorbed within the constant term. If 
the main halo accretes a subhalo of mass SM <C M 
during the time interval St <§C t, then the entropy 
of the entire system, including the pre-accretion 
entropy of the subhalo's baryons, rises by 

SS = S(M + SM, t + St)- S(M, t) 
- S(SM,t) 

fbSM 
jirrip 



M St , / M \ 
1 + TSM +ln {sM) 



(28) 



The first two terms in the lower expression corre- 
spond to the heat input needed to maintain self- 
similarity in the merged system, which presum- 
ably comes from thermalization of incoming ki- 
netic energy. The third term corresponds to the 
entropy associated with raising the subhalo's gas 
to the temperature of the main halo. The most 
efficient way to increase the entropy of that gas 
is to raise its temperature while keeping its den- 
sity constant, perhaps through conduction of heat 
from the main cluster facilitated by rapid mixing. 8 
However, the subhalo's baryons are more likely 
to be compressed or shock heated before mixing 
with the main intracluster medium, in which case 
a portion of the third term corresponds to ther- 
malization of incoming kinetic energy within the 
subhalo's gas. 

Smooth accretion achieves self-similarity in a 
different way. From equation (10) it is evident 



'Heat input that raises the temperature of a monatomic 
ideal gas from T\ to T2 at constant density increases its 
specific entropy s by (3/2) ln(T2/Ti). The temperature dif- 
ference between self-similar halos is T2/T1 = (M/8M) 2 / 3 , 
corresponding to a specific entropy increase ln(M/SM). 
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that cold accretion can produce self-similar clus- 
ters as long as ( = din M/ dint is the same 
for all clusters; the entropy profile produced by 
smooth accretion for constant £ is K(M g ,t) = 
K sm (t)[M g /f b M(t)] 2 / 3+2 / 3 i. In that case, the 
classical thermodynamic entropy of a cluster's 
baryons is 



S : 



fbM 
\im v 



In K 3 ' 2 

sm 



SO 



(29) 



which explicitly includes the constant term aris- 
ing from the cluster's internal structure. Smooth 
accretion of gas mass fbSM with pre-accretion en- 
tropy K\ <C K sra during a time interval St there- 
fore adds 



SS = 



fbSM 
fim p 



In 



(30) 



to the thermodynamic entropy of the whole accret- 
ing system, after the initial entropy of the accret- 
ing gas has been subtracted. The entropy increase 
in this case is due entirely to shock heating of in- 
coming gas at the accretion shock, and there are 
no additional heating terms. 

5.3.2. Implications of Accretion Heating 

The preceding analysis suggests that hierar- 
chical merging produces entropy exceeding the 
amount needed to raise the accreting gas to the 
temperature of the original halo. A detailed ac- 
counting of how this heating occurs is beyond the 
scope of this paper, but we would like to suggest 
a qualitative picture that can be tested and re- 
fined with numerical simulations. Here we con- 
struct a toy model for distributed heating owing 
to accretion and show how it affects the evolu- 
tion of the intracluster entropy distribution. If 
distributed heating is the dominant mode of en- 
tropy growth, then the entropy of gas deep within 
a cluster should rise with time, possibly compen- 
sating for some of the radiative cooling of the clus- 
ter core. 

Let the rate at which accretion deposits heat en- 
ergy into the intracluster medium be eT^fbM / \im p . 
If this energy deposition rate were equal to the 
rate at which kinetic energy passes through the 
virial radius, then we would have e = 1. However, 
accreting gas that is denser than gas at the virial 
radius plunges to smaller radii, gaining additional 



energy to be thermalized as it descends into the 
cluster. A dense gas blob eventually comes to rest 
when its entropy equals that of the surrounding 
gas. Thus, the total rate of heat input owing to 
lumpy accretion might be somewhat larger than 
the rate at which kinetic energy flows through the 
virial radius, allowing e to exceed unity. 

Distributing this heat input equally among all 
the cluster's gas particles leads to the following 
expression for the evolution of entropy with time 
within a Lagrangian gas parcel, including radiative 
cooling: 



d \ nK 3/2 ^ eT<j> d \ nM i 3 Kt' (T) 1 
It ~~ ~T~ dint 1 ~ 2 K 3 / 2 T 



■ (31) 



A more general model could be constructed by al- 
lowing e to depend on K, thereby accounting for 
inhomogeneitics in heat input. The mean value of 
e, integrated over time, and its spatial variations 
within a cluster would be interesting to measure 
in numerical simulations of hierarchical merging. 

For purposes of illustration, let us consider 
an extreme case of distributed heating explicitly 
tuned to preserve self-similarity. If T w and 
K c <C K, the cluster will remain self-similar as 
long as 

/dlnMV 1 

^ 1+ Unr) • < 32 > 

Figure 15 shows how the value of K associated 
with a Lagrangian gas parcel changes with time in 
a cluster of mass 10 14 Mq at to according to 
equation (31) with e set by equation (32), assum- 
ing the accretion history from Tozzi & Norman 
(2001). For simplicity, we have assumed T = 
throughout the cluster. The dotted lines indicate 
how K(t) of a gas parcel evolves after it appears 
in a certain position in the K-t plane. Gas pare- 
cels that have become incorporated into the clus- 
ter should evolve along these paths; however, the 
starting points of the paths are arbitrary because 
this model does not specify the entropy of a gas 
parcel just after it has accreted. The solid line 
shows K<f,(t), which runs approximately parallel 
to the neighboring K(t) tracks because we have 
chosen e to reproduce self-similar entropy growth. 

In this scenario, distributed heating owing to 
accretion is quite important early in the cluster's 
history. The long-dashed line labeled K t h corre- 
sponds to the K(t) track that ends at K(t ) = 0. 
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Fig. 15. — Change in entropy with time in the intracluster 
medium of the main halo owing to distributed heating and 
radiative cooling. The dotted lines show how Lagrangian 
gas parcels move through the K-t plane, according to equa- 
tion 31 with e = l+(d In M/d hit)- 1 and T = T^. The solid 
line indicates the characteristic entropy K$ (t ) of the grow- 
ing dark- matter halo. The short-dashed line gives K Bul (t) 
for £ = 0.5 and no preheating. The long-dashed line shows 
the threshold K t ^ below which a gas parcel is destined to 
condense by t = to. The dot-dashed line shows the cool- 
ing threshold K c (T,f,) corresponding to the characteristic 
temperature T^lt) of the halo. Because of the heat in- 
put generated by accretion, the entropy of the intracluster 
medium can increase with time, amplifying the effects of 
early entropy input and reducing the amount of gas that 
condenses by t = to. 

Gas blobs below this line are destined to cool and 
condense by the present time, unless internal feed- 
back intervenes. Notice that a relatively small 
amount of initial entropy (< 10 33 erg cm 2 g~ 5 / 3 ) 
can prevent intracluster gas from condensing by 
the present day, because distributed heating acts 
as an entropy amplifier. 

The situation is quite different in Figure 16, in 
which we have computed K{t) with e = 0. All the 
gas with initial entropy < 2 x 10 33 erg cm 2 g~ 5 / 3 
at early times condenses by the present time. If 
smooth, cold accretion is the dominant mode, then 
gas parcels enter the K-t plane on the short- 
dashed line labeled K sia , and everything that ac- 
cretes before 0.2io will radiate away all of its en- 
tropy. Substantial amounts of preheating or inter- 
nal feedback (> 2 x 10 33 erg cm 2 g~ 5 / 3 ) are neces- 
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Fig. 16. — Change in entropy with time in the intracluster 
medium of the main halo owing to radiative cooling only. 
The lines have the same meanings as in Figure 15. In this 
scenario, gas that accretes early in time is much more likely 
to condense by t = to- 

sary to prevent this gas from condensing. 

In both scenarios, the fate of a gas parcel ac- 
creting onto the cluster depends on where it enters 
the K-t plane as it becomes incorporated into the 
intracluster medium. For example, if the core of 
an accreting subhalo does not gain enough entropy 
during the merger to exceed K t h(t), then it will 
ultimately sink toward the center of the cluster 
and condense, unless some other heat source pre- 
vents it from doing so. The level of this threshold 
depends critically on e, which is why this param- 
eter should be measured using numerical simula- 
tions. Sampling of the core gas in the simulations 
of Navarro et al. (1995) shows that K (t) indeed 
rises with time (see their Figure 8), implying e > 0, 
but the data presented are insufficient for measur- 
ing the heating rate. 

If e is close to the similarity-preserving value as- 
sumed in equation (32) , then the origin of "cooling 
flow" gas in hierarchical merging differs substan- 
tially from the smooth accretion case, in which the 
lowest-entropy gas is the earliest gas to accrete. 
Large values of e strongly boost the entropy of gas 
accreted at early times, preventing it from cool- 
ing later on. Instead, K t h(t) reaches its maximum 
value near to/2, suggesting that the gas most likely 
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to be part of a cooling flow is the gas that accretes 
with low entropy at intermediate times. 

Unfortunately, the simple model we have devel- 
oped provides little insight into where a parcel of 
accreting gas should enter the K-t plane after a 
merger. Shocks probably raise the entropy of this 
gas somewhat before it mixes with the intraclus- 
ter medium, but the magnitude of this entropy 
increase is uncertain. Tracking the entropy of 
gas parcels during a numerically simulated merger 
would help to quantify this important parameter. 

6. Summary 

The observable X-ray properties of a relaxed 
cluster of galaxies arc entirely determined by the 
entropy distribution of its intracluster medium 
and the shape of the potential well that confines 
that medium. Because intracluster entropy is so 
fundamental, we have sought to understand its 
origin in terms of some simple analytical models. 
This analysis leads us to conclude the following: 

1. A classical smooth accretion shock produces 
an intracluster entropy distribution with a 
form much like the self-similar entropy pro- 
files observed in both massive clusters and 
in simulations of hierarchical structure for- 
mation. However, the normalization of the 
entropy profile arising from smooth accre- 
tion is ~2-3 times higher, depending on the 
accretion rate. Smooth, steady, isotropic 
accretion produces the maximum amount 
of entropy for a given M because smooth- 
ing of the incoming gas minimizes the mean 
mass- weighted density entering the accretion 
shock (§ 2). 

2. Preheating smoothly accreting gas to an en- 
tropy level K\ boosts the postshock entropy 
profile by an additive term w 0.84ffi, as long 
as the entropy of preheating does not exceed 
the characteristic entropy of the accreting 
halo (§ 3.1). However, preheating that ex- 
ceeds the characteristic entropy of the ac- 
creting halo strongly suppresses entropy pro- 
duction at the accretion shock because it in- 
flates the intracluster medium, pushing the 
accretion shock well outside the virial radius 
and thereby reducing the shock velocity at 
the accretion front (§ 3.2). 



3. Because smooth accretion of cold or mod- 
erately preheated gas produces ^2-3 times 
more entropy than observed in both massive 
clusters and simulations, inhomogeneity of 
the incoming gas must play a role in setting 
the entropy level of the intracluster medium. 
This difference in entropy production be- 
tween smooth accretion and hierarchical ac- 
cretion allows for an interesting mode of sim- 
ilarity breaking. If preheating has been able 
to eject gas from the small subhalos accret- 
ing onto groups, then accretion of baryons 
onto groups may be smoother than accretion 
onto clusters, enhancing entropy production 
at the accretion shock (§ 3.5). 

4. Smoothing of gas associated with low-mass 
halos might explain some interesting dif- 
ferences between the entropy profiles of 
groups and clusters. Isentropic groups can- 
not satisfy the observed L-T relation while 
having the observed core entropy ~ 100- 
150keVcm 2 (§ 4.1). Polytropic models of 
groups that satisfy both observational con- 
straints have entropy gradients implying en- 
tropy levels at <~ r 2 oo that substantially ex- 
ceed what hierarchical accretion can produce 
(§ 4.2). Direct observations indicate that 
entropy at the outskirts of groups is indeed 
^2-3 times higher than values derived from 
self-similar scaling of clusters (Finoguenov 
et al. 2002), perhaps owing to a transition 
from lumpy accretion to smooth accretion 
occuring on group scales (§ 4.3). 

5. Because the normalization of intracluster en- 
tropy appears to be sensitive to the density 
distribution of incoming gas, we would like 
to understand what sets that normalization 
in the case of hierarchical accretion. How- 
ever, a naive model in which all the incoming 
energy is thermalized by shocks within sub- 
halos having an average density similar to 
that of the main halo fails to produce enough 
entropy (§ 5.1). 

6. Given the failure of this simple model to ex- 
plain the observed normalization of the intr- 
acluster entropy distribution, we are driven 
to consider more complex scenarios for en- 
tropy generation. One possibility is that 
some of the accreting gas is not contained 



within subhalos, so that the incoming den- 
sity distribution is intermediate between 
smooth accretion and the assumptions of 
§ 5.1. Another possibility is that some 
of the incoming kinetic energy is thcrmal- 
ized within the existing intracluster medium 
through shocks and turbulence stimulated 
as accreting gas lumps circulate through the 
cluster (§ 5.3). 

7. If the kinetic energy of accretion is indeed 
an important source of distributed heating 
within a cluster, then the development of 
entropy in the lumpy accretion case is quali- 
tatively different from the smooth accretion 
case in which all entropy is generated at 
the accretion shock and is deposited exclu- 
sively into the accreting gas. One can test 
this possibility by measuring the rate of en- 
tropy growth deep within the accretion radii 
of simulated clusters (§ 5.3). 
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APPENDIX A 



SELF-CONSISTENT SHOCK RADIUS 



The position of the shock radius in an intracluster medium bounded by accretion pressure depends on 
both the mass accretion rate and the internal structure of the cluster. In this appendix, we develop a useful 
approximate solution for the shock radius in terms of the shock- velocity parameter £, based on a polytropic 
model for structure of the intracluster medium. The cluster potential is assumed to be of the NFW type, so 
that 

2T A ln(l + c A r/r A ) r A 



<Kr) = - 



(Al) 



pm p ln(l + c A ) — c A (l + c A ) 1 r 

where T A = GM A pm p /2r A and M A is the total mass within r A = r ta /2. We also assume that the pressure 
and density of gas in the equilibrium configuration are related by a polytropic equation of state, with 
P{r) oc [p(r)] 7off . Simulations, observations, and the smooth-accretion analysis of § 3.2 all suggest that 
7 c ff w 1.2 is an adequate description of the equilibrium state (e.g., Markevitch et al. 1998, 1999; Ettori 
& Fabian 1999; Finoguenov et al. 2001; Molendi & De Grandi 2002); however, it is important to realize 
that this relationship is not the actual equation of state of the intracluster gas. It is merely a convenient 
description of the global equilibrium state which ultimately depends on the entropy distribution K(M g ) and 
the shape of the potential well. 

Given these assumptions, the equation of hydrostatic equilibrium reduces to 



dT 

— = -pm p 



7eff 



1 



7eff 



dr 



which yields the solution 



T(x) 
p{x) 

P{x) 

9{x) 
9o(x) 

F(c) 



T A g(x) 

TAp 9 [g(a;)pff/(7eff- 1) 

pm p 

9o(x) +91 
2(7cfl- - 1) 



F(ca) 



ln(l + c A x) 



x = 



7eff ' x\n(l + c) 

ln(l + c) 
ln(l + c) -cil + c)- 1 
r/rA ■ 



(A2) 

(A3) 
(A4) 

(A5) 

(A6) 

(A7) 

(A8) 
(A9) 



The two constants of integration, gi and p g , are determined by the two boundary conditions of the self- 
consistent solution. First, the total gas mass within the dimensionless shock radius x ac = r ac /rA must 
equal the total mass of accreted gas i&M A . Second, the ram pressure of accreting gas at x ac must equal 
P{x). Furthermore, the entropy just within the shock radius must be consistent with the assumed polytropic 
relation. 

In the case of cold accretion, entropy consistency and the pressure-balance boundary condition are satisfied 
if T(x ac ) = pm p v ac /3 and p(x ac ) = Api. The temperature condition leads to 



9i 



|(2-, ac )-^L^ F(cA) Ml±£^ 

3 7eff ln(l + c) 



(A10) 



For typical cluster parameters of c A m 8 and 7 c ff ss 1.2, we find ffo(l) ~ 0.56 and gi ~ 0.11 when r ac = r A . 
Furthermore, we find gi = when r ac w l.llr A . Because go(x) declines with radius, we always have go 3> gi 
for x < x ac ] therefore, we elect to construct an approximate solution by setting g\ — 0. 
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To apply the gas-mass boundary condition, we take advantage of the density condition at x ac . Because 
we have set gi = 0, we can then write the gas-mass boundary condition as 



4/91 



= 3 



9o(r) 



ffoOac) 



l/(7eff-l) 



x dx 



Defining c a 



CA^ac/^A and recalling that £ = 1 — x ac /2, we therefore arrive at 

(i - W 3 



for the polytropic model, with 



and 



7( 7 



£poly 



Jo 



[6p/(7cff , C a 



Ml + cy) 
y\n(l + c) 

1/2 



1/(7-1) 



y dy 



3 A' m 



^dlnM 
dint 



(All) 



(A12) 



(A13) 



(A14) 



so that 4/»i = (/ 6 Ap cr )(20 



-V2,-3/2 = 



p. Numerical integration shows that the expression 7(1.2, c) w (c/8) 



0.45 is an excellent approximation in the interesting range 4 < c < 15. Note that the dependence of £ on 
the accretion rate is encapsulated in p and the dependence on internal cluster structure is encapsulated in 
I{l,c). 

The situation changes if the intracluster medium is iscntropic with an entropy level K\ that is not deter- 
mined by the shock at the boundary. For example, when preheated gas with K\ > 7T sm accretes onto the 
cluster through a weak shock, the entropy internal to the cluster is determined primarily by the amount of 
preheating. In the isentropic case, we can set 7 e ff = 5/3, yielding an internal pressure 



P(x) 



T A f b Ap c 
jimp 



El 



-3/2 



(A15) 



where Ka = TA(/um p ) 1 (fbAp a ) 2 / 3 . At the bounding radius x ac this pressure must equal the accretion 
pressure, implying that 



(2 X a c) 



1/5 T -1 



(A16) 



We can therefore write 



where 



9(x) = 



4F(c A ) 



ln(l + c A x) ln(l + CAX ac ) 
.Tln(l + c A ) x &c ln(l + ca) 

5 



Q 



Q 



4F(c A ) 



5(^ac) • 



(A17) 



(A18) 



In order to apply the gas-mass boundary condition, we need to integrate g 3 ^ 2 x 2 dx. We can approximate 
this integral by noting that ln(l + cx)/xln{\ + c) w x -2 / 3 at x > 0.3 to within 10% for c ~ 10. Because 
the outer regions of an isentropic cluster dominate the overall mass, the gas-mass boundary condition then 
becomes 

3/2 

xl c J(Q) , (A19) 





a 3 


"4F(c A )" 






5 
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where 



j(Q) = /V 2/3 -i + Q) 3/ V<% • 

Jo 



(A20) 



In the limit Q < 1, corresponding to a zero-pressure boundary, the integral can be done analytically to 
give J(Q) = [3r(3)r(5/2)]/[2r(ll/2)] = 0.0762. In the limit Q > 1, corresponding to a constant-pressure 
interior, the integral simplifies to iQ 3 / 2 . However, the most relevant range for our purposes is Q ~ 1, because 
we are interested in clusters whose shock-generated entropy is comparable to that produced by preheating. 
For these intermediate values of Q, numerical integration shows that J(Q) ~ Q/2 is a decent approximation. 
Applying this approximation produces 



(1 - &en)£ 



3/25 
isen;<,i S en 



-28/25 



3/5 r 



-i 3/10 



4F(c A ) 



-6/25 



Ki_ 



27/50 



(A21) 



As an accreting cluster makes the transition from the isentropic to the shock-dominated regime, we need 
a hybrid solution. Let be the fraction of accreted gas that is isentropic at time t, let Xk be the 

dimcnsionless radius within which the gas is isentropic, and define §k = 9(xk)- The gas-mass condition for 
the polytropic (7 e gf = 1.2) region at x > Xk is then 



(1 - fx) 



9 



x dx 



where g a 



4 Pi J Jx K V5a 

g{x SLC ). The corresponding gas-mass condition for the isentropic region is 



(A22) 



/ 



K 



= 3{ 9K 



-) 

9k J 



3/2 



x dx 



4 Pl J \5ac/ Jo 

Combining these expressions yields an expression for £ analogous to that for the pure polytropic case: 



(A23) 



1 



(6p)2 



9 

5ac 



y dy 



9k 

.9ac 



\ 3/2 

JL] 



9k J 



y dy 



(A24) 



with y = x/ 'x ac . The first integral approaches 7(1.2, c) as 
integral, we note that 9k ~ 9o{xk) at x > Xk implies that 



0, as long as g w g - To estimate the second 



12 



-9k 



\x K ) 



-2/3 



- 1 + 



12 



for x < xk- Thus, the second integral approaches w (12/5) 3 / 2 J(5/12) as 
following approximation to compute ^ for the transitional case: 



(A25) 



x ac . We therefore use the 



(1 



1 



(6/5) 



^ \l - f K ) ■ 7(1.2, c) + / x • (12/5) 3 / 2 J(5/12)l 



(A26) 



The solid lines in Figure 7 show ^ = mm(£ po i y , ^tran S) 6sen) determined using the same cluster parameters 
that Tozzi & Norman (2001) used to compute entropy models for clusters of 10 14 h~ x M Q and 10 15 h~ x M Q . 
We converted their dM/dz relations for ACDM to dM/dt relations, yielding M cx f0.93i-i.oiog/r < f or ^ 
IO^Ii^Mq cluster and M cx £0.547-0.914 log if t for the lQ^h~ l M & cluster. In order to reproduce the 
time-dependence of the concentration parameter we interpolated a power law in (1 + z) between their mass- 
concentration relations at z = and z = 1, giving c A = 8.5(M A /10 15 hr 1 M Q )-°- 086 (l + z)-°- 65 . We 
also neglected the difference between c ac and c A in computing I(j,c). Note that this model is not a good 
representation when £ <C 0.5 because it assumes that the time when gas reaches the shock front is twice the 
time it took to reach its turnaround radius. 
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APPENDIX B 



INHOMOGENEOUS ACCRETION AND ENTROPY PRODUCTION 

The following computation explicitly demonstrates that inhomogenous accretion produces less entropy 
than homogeneous accretion. Define a dimensionless baryon density x = Airr 2 vp/M and let f(x)dx be the 
fraction of the accreting volume with density between x and x + dx, so that J xf(x)dx = 1. The fraction 
of accreting mass with density in this range is then xf(x)dx, because of how density is defined. Thus, the 
mean mass-weighted density exceeds the mean volume-weighted density by a factor 

J x 2 f(x)dx = J(x- l) 2 f(x)dx + 1 . (Bl) 

In the homogeneous case, f(x) is a delta function at x — 1, and the mass- weighted density equals the 
volume-weighted density. However, if there is any inhomogeneity, the integral on the right-hand side will be 
greater than zero, implying that the mass-weighted mean density must be larger than the volume-weighted 
mean density. 

We are interested in the mean mass-weighted entropy of gas accreting at time t. If all the accreting gas 
moves with the same velocity and passes through an accretion shock at the same radius, then the mean 
mass-weighted entropy of accreting gas is 

K(t) = K sm {t) J x- 2 / 3 • xf(x)dx = K sm (t) J x 1 /\f(x)dx . (B2) 

We define y = x 1 / 3 and g{y)dy = f(x)dx, so that / y 3 g(y)dy = J g{y)dy = 1, implying / (y 3 - l)g{y)dy = 0. 
This equation leads to 

J(y- l)g(y)dy = -lJ( v + 2)(y- lfg{y)dy . (B3) 

Because the integrand on the right-hand side is always positive, the integral on the left-hand side must be 
less than zero, except in the homogeneous case. Thus, we have 



x l ' s f(x)dx= / yg(y)dy<l , (B4) 



K(t) 
K sm {t) 

demonstrating that homogeneous accretion maximizes entropy production in an accretion shock. Notice that 
this conclusion also applies to smoothing of the flow in solid angle or in time, meaning that the flow must 
also be isotropic and steady in order to maximize the postshock entropy. 
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